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Abstract 

State preparation via conditional output measurement on a beam splitter is stud- 
ied, assuming the signal mode is mixed with a mode prepared in a Fock state and 
photon numbers are measured in one of the output channels. It is shown that the 
mode in the other output channel is prepared in either a photon-subtracted or 
a photon-added Jacobi polynomial state, depending upon the difference between 
the number of photons in the input Fock state and the number of photons in the 
output Fock state onto which it is projected. The properties of the conditional 
output states are studied for coherent and squeezed input states, and the proba- 
bilities of generating the states are calculated. Relations to other states, such as 
near-photon-number states and squeezed-state-excitations, are given and propos- 
als are made for generating them by combining the scheme with others. Finally, 
effects of realistic photocounting and Fock-state preparation are discussed. 



1 Introduction 



Over the last years numerous workers have studied various nonclassical states 
of radiation and proposed schemes for producing them. Particular interest has 
been devoted, e.g., to Fock states (for a review, see [[J) and states derived from 
Fock states by coherently displacing and/or squeezing them |2], [| [|, ||, |6|, super- 
positions of mesoscopically distinguishable states, such as Schrodinger-cat-like 
states (for a review, see [?j), near-photon- number states (also called crescent 
states) || ||, [10], [11], 0, [0|, binomial states |14], [15|, inverse binomial states 
[P"7| , [TBI , squeezed-state excitations [|Tj], |2D| and the SU(2) and SU(1,1) minimum- 



uncertainty states |T0| , |21| , |22| . Another interesting class of nonclassical states 
that have been a subject of increasing interest are photon-added and photon- 
subtracted states that are obtained by repeated application of photon creation or 
destruction operators, respectively, on a given state 



29, 30 



Similarly, states obtained by the repeated application of the inverse boson oper- 
ators have also been considered |3l[ . 

Despite the large body of work, only few of the above mentioned nonclas- 
sical states have been generated experimentally so far. Designing of realistic 
schemes for generating specific quantum states and realization of the schemes in 
the laboratory have been one of the most exciting challenges to the researchers. 
A promising method of quantum state engineering has been conditional mea- 
surement, e.g., generation of a desired state by state reduction in one of two 
entangled quantum objects owing to an appropriate measurement on the other 
object. Typical examples that have been considered for nonclassical state genera- 
tion via conditional measurements are the interfering fields in the output channels 
of a beam splitter |23|, ^6], [3Sj, |33|, waves produced by parametric amplifiers 

~~8j, |39|, [|(J and degenerate four-wave mixers |U| fill 

or 
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pi m, m 



and systems of the Jaynes-Cummings type in cavity QED [| 
trapped- ion studies ]4(| [47], Further, state reduction via continuous mea- 
surement has also been considered ED], ISO, FTD, B31 EM. 



In this paper we study the class of states generated by conditional output 
measurement on a beam splitter in the case when an input mode prepared in 
some quantum state and another input mode prepared in an n photon Fock state 
are mixed and in one of the output channels of the beam splitter a photon-number 
measurement yields m photons. We show that the conditional output states are 
photon-subtracted (n < m) or photon-added (n > m) Jacobi polynomial states, 
i.e., states that are obtained by (\n — m| times) repeated application of either 
the photon destruction operator or the photon creation operator, respectively, 
to Jacobi polynomial states. It is worth noting that the scheme can be used 
to generate photon-subtracted and photon-added Jacobi polynomial states for 
various classes of input quantum states, such as thermal states, coherent states, 
squeezed states and displaced photon-number states. In particular, for n = 
and m = 0, respectively, ordinary photon-subtracted and photon-added states 
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H m, H, |2|, H H H, are observed. 

In order to illustrate the nonclassical properties of photon-subtracted and 
photon-added Jacobi polynomial states, we study them for coherent input states 
in more detail. We analyze the states in terms of the photon-number and 
quadrature-component distributions and the Wigner and Husimi functions, and 
we calculate the probability of producing them. Further, we briefly address the 
case of squeezed vacuum input. It is worth noting that in this case the produced 
photon-subtracted and photon-added Jacobi polynomial states - similarly to or- 
dinary photon-subtracted and photon-added squeezed vacuum states [^, |25|] 
are examples of classes of Schrodinger-cat-like states. 

We further study the relation of the conditional output states to other classes 
of nonclassical states. In particular we show that near-photon-number states 
[[| 0, |TD|, [II], [12], |TB[ and squeezed-state excitations JT9], |2U| can be generated by 
photon adding and subsequent coherent displacement and/or squeezing. It also 
turns out that photon-subtracted and photon-added Jacobi polynomial coherent 
states are finite superpositions of ordinary photon-added coherent states, which 
for themselves are finite superpositions of displaced Fock states ||28|| . Similarly, 
photon-subtracted and photon-added squeezed vacuum states can be regarded as 
two different finite superpositions of squeezed number states. 

This paper is organized as follows. Section presents the basic scheme for gen- 
eration of photon-subtracted and photon-added Jacobi polynomial states. The 
properties of the states are studied in Sections [3] [5], with special emphasis on 
coherent input states (Section |J) and squeezed vacuum input states (Section ||). 
Relations to other states are given in Section ||. In Section [7| effects of nonperfect 
preparation and measurement of photon-number states are addressed. Finally, a 
summary and concluding remarks are given in Section || 



2 Scheme of conditional measurement 

Splitting and mixing optical fields on beam splitters are basic manipulations in 
classical as well as in quantum optics. The input-output relations of a lossless 
beam splitter are well known to obey the SU(2) Lie algebra |55|, |56|. In the 



Heisenberg picture, the photon destruction operators of the outgoing modes, b k 
(k = 1,2), are obtained from those of the incoming modes, as 

2 

bk = Tk,k' Ofe'j (1) 
fc'=i 



where 

(T, ,.\ — J<po ( 

— sin9e llpR cos 9 e 



{J-k,k')—d | _ -„ a a -i VR „ na Q a -upr \ z ) 
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is a SU(2) matrix whose elements are given by the complex transmittance T and 
reflectance R of the beam splitter, 

T = cos # e iLpT , R=sm6 e iLpR . (3) 

In the Schrodinger picture, the density operator is unitarily transformed, whereas 
the photonic operators are left unchanged. In this case the output-state density 
operator g OVLt can be related to the input-state density operator g in as 

£out = V^QmV, (4) 

where V can be given by ]55], [56 1 



y — g-iCvT-Vii)^ e -2i8L 2 e -i(v>T+V>R)L3 ^ 

with 

^2 = Ti(o\a 2 - a\ax), L 3 = \(a\ai - a\a 2 ). (6) 

Note that <po is a global phase factor, which may be omitted without loss of 
generality, <po = 0. Applying elementary parameter-differentiation techniques 
0, we can derive the operator identity 

— 2i9L2 tanSa^ai 2 In cos L3 — tan a j S2 fy^ 

which [together with equation (|3|)] enables us to rewrite V', equation (|5|), as 

•y\ _ ijfhi g— -R*o|oi gRaJa2 rp~n 2 ^g^ 

where n fc = a|.afc. 

An outline of the experimental setup is depicted in Figure III. A field mode 
prepared in a state described by the density operator g inl is mixed at a beam 
splitter with another mode prepared in a Fock state \n). The input-state density 
operator can then be written as 

kn{n) = ftni ® \n) 2 %{n\. (9) 

Using equations @, (§) and @, after some algebra the output-state density 
operator g OVLt = g out (n) can be given by 



oo oo n n 



l J I Z=0 m=0A;=0 j=0 



n\(n\{ n — k+m\( n—j+l\ 

vmm\\k)[j){ m )[ i ) 

xT" 1 a™(a t 1 )^ inl a{(a t 1 ) / (T*)" 1 ®|n-A;+m) 22 (n-j+/|. 



(10) 
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From equation (jlOf) we see that the output modes are highly correlated to each 
other in general. When the photon number of the mode in the second output 
channel is measured and m photons are detected, then the mode in the first output 
channel is prepared in a quantum state whose density operator £ outl (n, m) reads 
as 

* , n n _ 2{m\Q ont {n)\m) 2 ^ 
° Utl ' Tri{ 2 (m|£ ou t(n)|m) 2 } ' 

The probability of such an event is given by 
P(n,m) = Tri{ 2 (rn\g ont (n)\m) 2 } 

lRl ^ nl ff(-\R\ 2 ) u+k) ( m W m ^ 



\T\ 2m m\ , \j — vj\k — uj 

" n\\T\ 2n fn+j\ (n + k\. .... 

where the abbreviations 

v = n — m, fi = max(0, z/), 5 = fi — v (13) 

have been used. Let us now assume that the mode in the first input channel is 
prepared in a mixed state 

ftnl = E^I $ )( $ l ( 14 ) 

(J2<s>P<s> = 1, <p$ < 1). Combining equations ( |T0D and ([□]) and using equation 
(|14D , we find that the mode in the first output channel is prepared in a state 



Q out l{n,m) = J2P$\ ^n,m)(^n,m\, (15) 

where 

•N"n,m being the normalization constant, 

_ » » (-\R\ 2 ) k +i (n\ (n\ 



X 



($|(T*) ftl a|(ai) J '- l/ a^(a t 1 ) fe r ftl |$). (17) 
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3 Photon-subtracted and photon-added Jacobi 
polynomial states 

The properties of the conditional output state essentially depend on whether 
photons are effectively subtracted (n < m) or added (n > m) . To be more specific, 
from equation flltf ) we obtain for n < m the output state 

\_ a/--i/2 (~\ R \ 2 ) k ( n ia^ia fc rgti fc T"i^ 



fc= 



where the notation 

a fc (a t ) fc = (n + l)(n + 2) • • • (n + k) = k\ (19) 

has been introduced (n = hi). For n>mwe derive [on using the relation n(a^) 1 
= (tf) l (n + l)] 



|*n,m> 



2\i/ 



m!A/"n,m 
m f I D|2\fc 



'-|i?| 2 )^! 



mWr 

X 



1/2 



m / 



From equations flI8|) and ( |2"0] ) it can be shown (Appendix A) that the conditional 
output state is of the form 



I * n,m/ 



' aWp(H«-™)(2|T| 2 -l)T™|$) forz/<0, 

(21) 

(at) t/ P^"- m )(2|T| 2 -l)T™|<l>) fori/>0, 



where P| a ' /3 ' ) (z) is the Jacobi polynomial. The following procedure is seen to 
yield the conditional output states from a chosen input state |$) = J2kLo c k\k). 
(i) Replace the Fock-expansion coefficients with c' k ~PS 1 fe ~ (2|T| 2 — 1) T k Ck 
to obtain a state |$'). (ii) Subtract photons from the state |<E>') by repeated 
application of the photon-destruction operator to it or add photons to state |$') 
by repeated application of the photon-creation operator to it. In what follows 
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we will refer to the states |$') as Jacobi polynomial (JP) states (in analogy to 
Hermite polynomial and Laguerre polynomial states |5^|). Note that for 
typical classes of states the input state |<3>) and the state \ty) ~T n |<3>) belong to 
the same class of states P5| . We see that the conditional output states \^f n ,m) 
produced in the scheme can be regarded as photon-subtracted Jacobi polynomial 
(PSJP) states (n < m) and photon-added Jacobi polynomial (PAJP) states (n 
>m). It should be pointed out that the PSJP states and the PAJP states are 
essentially different from each other in general, because of [a, a*] 7^ 0. Clearly, 
when the input state is a Fock state, |<&) = \ k), then the conditional output states 
are the Fock states \k + n — m). Let us mention that when m = n (i.e., v = 0), 
then 

l*«,n> = K,W A - n) (2|T| 2 -l)T ft |$) . (22) 



4 PSJP and PAJP coherent states 

To treat the states in a unified way, let us return to equation ( ]T5| ) and first 
consider Glauber coherent input states 

00 ok 

m = \{3) = e-WY. J 7u\k), (23) 
with (3=\/3\e i(p i 3 . Equation ([TBD then reads 



^^ST^fC)^™' (24) 



where (3' = T(3 and Af^ m = ''Af n ,m- Applying standard operator ordering 
techniques p0| , we may write 

min{m,n} / \ | 

a m (a f r= E (7j(^oI (at) ^ am ^' (25) 

and hence 



• n.m k fi 

xg(^C)(^)'l^>, (26) 
from which the Fock-state expansion of \^ n ,m) can easily be obtained to be 

e -|/3'| 2 /2 n / \ 

y-' v n,m fc=/i \ / 

x E , 77^£^V^+^b + ( 27 ) 
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(for the photon- number statistics, see Appendix B). Using the identities |JI 



f^ T{l + u)\lJ T(n + v) 

and 

" g Lf(x) _ (l+t)"Lg[tx/(l+t)] 
^ r(a + Z + l) (ra-Z)! T(a+n+l) 1 J 



[L^(^) is the associated (or generalized) Laguerre polynomial we may give 
equation (^) in the more compact form of 

\T\ 2n n\ \-\R\ 2 /\T\ 2 T 

I ^ ra,m/ 



IT 



x 1^ (S^otWfiW). (30) 



In a similar way it follows that 



M' = — ^- TT(-\R\ 2 ) k+j ( n )( n ^ 



where 



(1) J fc!o*-*L£-*(-|a| 2 ) for/>fc, 

j " 1 l\(a*) k - l L?- l (-\a\ 2 ) for Z < jfe. 1 J 



From equations (|26|) or fl3~0|) we find that PSJP and PAJP coherent states are 
finite superpositions of photon-added coherent states. In particular for n < m 
equation ( PD| ) reduces to 

(33) 

from which we see that when \f}'\ Cl and 7 = (I-Rp/ITI 2 )/?' finite, then the PSJP 
state \^ n ,m) is (approximately) a superposition of n Fock states, because of \j3') 

«|0>- 

The probability of producing PSJP and PAJP coherent states can be obtained 
from equation (|12|). After some calculation we derive (Appendix C) 

\R\~ 2u n\ 



\T\ 2m m\ 



n n 

X 



k=fj,j=ii 
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where 



E 

1=0 



L;":'(-i«i 2 )ia 



2/ 



for 1/ > 0, 



E 

1=0 



(35) 



LJ+i (-|«| a )l«| 2M forz/<0. 



In Figure |2] examples of P(n, m) are plotted for two absolute values of the beam- 
splitter transmittance. We see that P(n, m) is an oscillating function of the 
absolute value of the coherent input amplitude, which is due to the interference 
of the incoming fields at the beam splitter. As expected, for n<m the probability 
P(n, m) goes to zero as the coherent amplitude does. 

The superposition of photon-added coherent states, equation (|26)) , gives rise 
to strong quantum interference, as it can be seen from the quadrature distribution 



p n (x,ip\m) = |(x,y?|^ nim )| 2 . 



Using the Fock-state expansions (27) and 

\x,<p) 



7T 



l/4 e -x 2 /2 e 



H fc (x)|fc) 



(36) 



(37) 



[H fe (:r), Hermite polynomial] of the states |^ n ,m) an d , respectively, and 
recalling the identity |JT] 



oo z k 



E TJ H fc+n(^) = exp [2xz - z 2 J R n (x - z) 



k=0 



(3? 



we find that 



p n {x, (p\m) 



r\-2u 



x exp< — 



ML 



x—y2\ft'\ cos((p — (p/3') 



EH* 

k=fi 



2\k 



n' 
. k , 



E 

i=fi 



'k\ (2-l/Te^)' 



Hi(x-2-^/3'*e^) 



(39) 



where <^g/ = <^g + Examples of p n (x, y?|m) are plotted in Figure |3| for ft' = 2.07 
and (a) n = 2, m = 3 [P(n, m) »s 9.7%] and (b) n = 3, m = 2 [P(n, m) « 6.7%]. 
From the figure it is clearly seen that the states are extremely non-Gaussian 
squeezed coherent states owing to quantum interference. 
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Next let us consider the Husimi function 



Q n (x,p\m) = — |(a|^ n , m )| 2 , 

Z7T 



(40) 



with | a) being a coherent state and a = 2 1 / 2 (x+ip). Expanding \a) and \^ n ,m) 
in the Fock basis, on applying equations (p3|) and p7|) respectively, we derive 



Q n (x,p\m) 



exp[— 


a- f3' 


2 ] 


2tt|/5' 


2u j\f' 

J n,m 



X 



k=fi 



2\k 



n 

. k , 



(41) 



We again use the relations (|28| ) and (^) and rewrite equation ([11]) as 



Q n (x,p\m) 



\T\ 4n (n\) 2 (-\RTal3'\) 2 » 



x e 



| a -/3'| 2 



L 



n + S)\ 

M 



n— /i I 



|i2 


2 a/T 




T 


2 



(42) 



Figure £| shows plots of the Husimi function for (3' = 2.07 and (a) n — 2, m = 
3 and (b) n = 3, m = 2. It should be pointed out that (the square root of) the 
difference between the height of the Husimi function and 1 / (2ir) can be regarded, 
in a sense, as a measure of the degree of nonclassicality |J2[. For the most 
classical states, i.e., for the Glauber coherent states and only for them, the Husimi 
function attains the maximum height of l/(27r) ~ 0.16. Thus from Figures ^(a) 
and ^(b) we see that the PSJP coherent state (n < m) is more classical than the 
PAJP coherent state (n > m) - a behaviour that is observed for a wide range of 
values of \(3'\. 

Among the phase-space functions that have been shown to be inferable from 
measurable data, the Wigner function 



W n (x,p\m) = - [ dye 2ipv {x-y\^ n ^){^ n)m \x+y) 

71 J 



(43) 



(\x) = \x,(p)\ 9= o) reflects quantum features most distinctly. In order to calculate 
it, we again use the Fock-state expansions ( p7| ) and ( p7|) [together with equation 
After some calculation we derive 



W n (x,p\m) = 



c 2 + |( / 3'4 7 3'*) 2 -v / 2x(/3'+/3'*) 



•i\2v 
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(2 ^ 2 (3') l (2 1 / 2 (3'*) 1 ' r 2 +y[2ip+V 2 i/3 ,*_ 0/)] 
(l-v)\(V -v)\ J dye 

— OO 

x H i (a;-y-2- 1 / 2 /3 / )H,(x + 2 /-2- 1 / 2 /3 / *). (44) 

We then use the integral identity [ plfl 

+00 



J dxe x B. m (x+a)H n (x+b) 



hm—n 1 

and perform the ^/-integration to obtain 



2ab) (m<n), (45) 



-|a;+ip-v/2/3'| 2 n n , 

W n (x,p\m) = E E {H*f) fc+J 

I /lyv n,m k=fij=tJ, K 



i=fi i'=fi 



l)\V) {I - u))\{V - u)\ 
x x [, 3 l[v / 2(x+^-2- 1 /V)]}, (46) 



where 



{3) f (-1) W^L^af) for I > fc, 

Xi,kW <y (-iyi\(a*) k - l L?-\\a\ 2 ) for / < k. 1 J 

In Figure [5] plots of the Wigner function W n (x,p\m) are shown for /3' = 2.07 and 
(a) n = 2, m = 3 and (b) n = 3, m = 2. We see that W n (a;,p|m) for n> m is more 
structurized owing to quantum interference and stronger negative than for n < 
m, which again reveals that the PAJP coherent state is more nonclassical than 
the PSJP coherent state. 



5 PSJP and PAJP squeezed vacuum states 

Let us briefly comment on PSJP and PAJP squeezed vacuum states that are real- 
ized when the input state is a squeezed vacuum ,S(£)|0), with S'(^)=exp{ — ^(a*) 2 
— £*Q 2 ]} being the squeeze operator. Hence we may write 

1 00 \(2kVW 2 

|$) ee <?(O|0> = (1 - M 2 ); E * k \2k), (48) 
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where the notation k = /«(£) = — e t{f z tanh |£| has been introduced (£ = |£|e 11 ^ 
According to equation^), the conditional output states are then given by 



|*n,m) 



E 



\R 



2\k 



n 



(49) 



with Af^ m 
reads as 



[(1 — |/t'| 2 )/(l — I^H^^A^m,. In the photon-number basis |*n,m) 



I ^ n,m/ 



ML 



E 

k=fi 



\R 



2\k 



n 



(k-u)\ \k. 



X 



(p-z/ + fc)!±[l + (-l)P-"] 



\p), 



(50) 



with n' = T 2 k. From equation ( j5"0|) we easily see that when the difference between 
the number n of photons in the second input channel of the beam splitter and the 
number m of photons detected in the second output channel, i.e., the parameter v 
= n — m, is even (odd), then the mode in the first output channel is prepared in a 
PSJP or PAJP squeezed vacuum state \^f njm ) that contains only Fock states with 
even (odd) numbers of photons. Similarly to ordinary photon-subtracted and 
photon-added squeezed vacuum states [|3|, |24|, it can be shown that the PSJP 
and PAJP squeezed vacuum states are Schrodinger-cat-like states. In particular, 
from equation (|50"D it can be found that \^ n ,m) can be given by a superposition 
of two mesoscopically distinguishable states, 



where 



with 



#(+)\ + |#(-)\ 
n,ml 1 I n,ml ) 

1 oo 

EcS>') Ip>, 



AC 



(±> 



(51) 
(52) 



2\k 



n 



(p — + &)! 



r 



§(p-i/)+i 



p—v 



(53) 



A detailed analysis can be given in a way similar to that in |23|, p4]]. We therefore 
renounce the somewhat lengthy calculations here. 
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6 Relations to other states 



From equation (21) and the property of the Jacobi polynomials that 



P<«.»(1) = (-1)<P<"'«(-1) = ('+ Q j (54) 

we see that for sufficiently small values of |T| (|T| — >• 0) or values of |T| close to 
unity (\T\ — > 1) equation (|2T|) approximately reduces to 

f a H T™|$) for v = n -m < 0, 

|*n,m) ~ I (55) 
[ (a)yT h \§) for v = n - m > 0. 

As expected, in these limiting cases the produced conditional states \^ n ,m) reduce 
to ordinary photon-subtracted [y < 0) or photon-added {y > 0) states. 



Further, equation (|2TD reveals that when the second input mode is in the vac- 
uum state and in the second output channels a nonvanishing number of photons 
is detected (^ = 0, m>0), then usual photon subtraction is observed, 

|*o,m> = (l/m!)AC /2 a m n*>, (56) 

independently of the value of T. Similarly, when a Fock state \n) is fed into the 
second input channel and a zero-photon conditional measurement is performed 
in the second output channel {n > 0, m = 0), then the mode in the first output 
channel is prepared in an ordinary photon-added state 

|* n ,o) = {-\R\ 2 ) n N^ /2 {cJ) n T h \<$>). (57) 

Photon-subtracted states and photon-added states of the type given in equations 
( |56| ) and (0), respectively, have been studied by several authors [^3], ^5], p6 
E7|,EaE9LE0|. 



6.1 Displaced Fock states 

Since the coherent states \/3) = D(j3)\0) [D{/3) =exp(f3a^ — f3*a)\ are the eigenstates 
to the destruction operator (a\/3) —/3\(3)), it is obvious that subtracting photons 
from a coherent state yields again a coherent state. Photon-added coherent states 
are highly nonclassical states. From 

(at)»|0) = (at)*£>G0)|O> 

we find that 

(at)»|/J) = D(j3){tf + F) n \0) 

= £(fyVfi.((3T- k D(m. (59) 
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Equation (p9|) reveals that photon-added coherent states are finite superpositions 
of displaced Fock states D(/3)\n) |2?| (for the properties of displaced Fock states, 
see 0, |3|, Hj. From equations (|26"D and (|3CiD we know that PSJP and PAJP coherent 



states can be given by finite superpositions of photon-added coherent states. 
Expressing the latter, according to equations (|58D and (^), in terms of displaced 
Fock states, we see that PSJP and PAJP states are also finite superpositions of 
displaced Fock states. 

Similarly, displaced Fock states can be given by finite superpositions of photon- 
added coherent states. To show this, we write 

D{(3)\n) = D(f3)^=(tfr\0) 



and hence 

D((3)\n) 



' b(P)(a^ n b\p)b(P)\0), (60) 




aj\f3). (61) 



Note that from equation ( ]6"T| ) and equation ( |3T)D for |T| 2 = \R\ 2 = 0.5 and (3 ~ 
1 it follows that |^n) ~ D(/3)\ n)\ n= \, i.e., the conditional measurement schemes 
realizes a coherently displaced single-photon Fock state. 



6.2 Near-photon-number eigenstates 

Near-photon-number eigenstates are an example of minimum uncertainty states 
that are defined by the eigenstates of the operator Y=n — i\f3\x{<p), which is asso- 
ciated with the "simultaneous" measurement of photon number and quadrature 
components |§. The states are also called crescent states and have been studied 
in a number of papers [§, |], |1(], [LJ], [12], [13| . In particular, they can be expressed 



in the form of 

M~ (a f + /FTO, (62) 

or alternatively, in terms of nonunitarily shifted Fock states, and it was shown 
that they can be generated by state reduction via photon-number conditional 
measurement in nondegenerate parametric down conversion ]13 . 
From equation ( |62"D it is easily seen that 

M ~ &V3)@) n D(PM = b(-(3)(a^ n \2(3). (63) 

Equation ( |6"3] ) reveals that near-photon-number eigenstates are coherently dis- 
placed photon-added coherent states, which offers the possibility of producing 



13 



them by conditional measurement on a beam splitter and subsequent coherent 
displacement. As it is depicted in Figure || a mode prepared in a coherent state 
| a) is first mixed with a mode prepared in a Fock state \n), and a zero-photon 
measurement is performed in one output channel of the beam splitter to prepare 
the mode in the other output channel in a photon-added coherent state (a^) n |2/3) 
(2/3 = T a). To realize the coherent displacement D(—(3), this mode and a strong 
local-oscillator mode prepared in state |ol) are then superimposed by another, 
unbalanced beam splitter of high transmittance T and low reflectance R such 
that (3 = -(R/f)a L . 

6.3 Squeezed Fock states 

Next let us consider a photon-added squeezed vacuum state, 

(fit)"£(0|o> = s(0&(0(tfrs(m 

= (l-\K\ 2 )- n/2 S(£)(a< + K*a) n \0) (64) 

[S(0 = exp{-I[£(at) 2 - Ca 2 }}, « = «(£) = -e^ tanh £ = |£|e^]. Using 
standard ordering techniques for boson operators we derive 



(at + ea)n = ^(l) (e)n "' 

[3] foAk 



x Hlfr( H ^ |i,rr ' (65) 



which enables us to rewrite equation ( 64]) as 

(atr£(0|o) = (i-NT" /2 

x E - , /, J = S(£)\n-2k). (66) 
fc =o 2 k k\^j(n-2k)\ 

Analogously, for the photon-subtracted squeezed vacuum states we find that 



a m 



5(0|0) = 5(0^(0^(010) 

= k(1- |K| 2 )- m/2 5(O(a t + ^ 1 a) m |0), (67) 



and hence 

a m 5(0|o) = K(i-|K| 2 )- m/2 



\ £ . 2>wl- a / (0|m " a> - <K 

fe=0 rh'.y yiii £1%)'. 
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Equations (|66"D and (|68[) show that photon-added and photon-subtracted squeezed 
vacuum states can be given by finite superpositions of squeezed Fock states 
S(^)\k), and it is worth noting that the two classes of states realize two classes 
of Schrodinger-cat states P3|, The extension of equations (^) and ( |B^ ) to 
photon-added squeezed coherent states (a^) n D(j3)S(£) |0) and photon-subtracted 
squeezed coherent states a m D(P)S(t;) |0) is straightforward. These two classes 
of states can be given by finite superpositions of displaced squeezed Fock states 
D(P)S(!;)\k) (for the properties of displaced squeezed Fock states, see [|, §). 
Displaced squeezed Fock states D(P)S(C,)\n) can be rewritten as 

D(P)S(0\n) = (nir^D^SiO^nO) 

= [(1 - H 2 )^!]- 1 / 2 ^/?)^ - K*a) n S{O\0). (69) 

From equations ( |69"D and (|65|) it is seen that displaced squeezed Fock states cannot 
be given by finite superpositions of the photon-added squeezed states in general. 
Note that the displaced Fock states D(P)\k) are finite superpositions of photon- 
added coherent states [see equation flBIp]. 



6.4 Squeezed-state excitations 

Squeezed-state excitations 

\(3,n;O~D(P)(tf + K*a) n S(O\0) (70) 

were introduced for diagonalizing the complete Gaussian class of phase-space 
functions [19, 20fl . Note that although equation (|70|) formally resembles equation 



squeezed-state excitations are quite different from displaced squeezed Fock 
states in general. It can be easily seen that 

~ D(^(o(^rm)\o) (7i) 

[cf. equation (0)]. We see that squeezed-state excitations are nothing but 
squeezed and subsequently displaced photon-added squeezed vacuum states, which 
implies the scheme in Figure |7| for producing them. Using a beam splitter, a mode 
prepared in a squeezed vacuum S'(C)lO) is first mixed with a mode prepared in a 
Fock state \n) and a zero-photon measurement is performed in one output chan- 
nel of the beam splitter to prepare the mode in the other output channel in a 
photon-added squeezed vacuum state (a^) n S(2^)\0) [k(2£) = T 2 k(()]. This state 
is then squeezed (with — £ as squeezing parameter), e.g., by degenerate para- 
metric amplification or by state reduction via quadrature-component conditional 



measurement in nondegenerate parametric down conversion ||34|| . Superimposing 
the outgoing mode and a strong local oscillator by an unbalanced beam splitter, 
the coherent displacement can be realized [j3 = (R/T)aL, cf. Section |0| . 
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7 Realistic experimental conditions 



Let us first address the problem of realistic photon detection. Unfortunately, 
there are no highly efficient and precisely discriminating photodetectors available 
at present. To overcome this difficulty, photon chopping ||§ was suggested for 
measuring the photon-number statistics. Let us remember that in such a scheme 
the mode to be detected is fed into an input channel of an optical 2iV-port array 
of beam splitters, the other N — 1 input ports being unused. Highly efficient 
avalanche photodiodes in the N output channels are used in order to record the 
coincidence event statistics. Since they only distinguish between photons being 
present or absent, the probability of obtaining k clicks when m photon are present 
is given by p3 

PNMm) = Y,PN(W)Mi, m tn) (72) 
i 

(rj is the detection efficiency), where 

P»Wm) = (*) £(-!)< ( k \k- l) m (73) 



for k<m, and P^(k\m) = for k> m. The matrix 



MUv)={^Jv l (l-v) m - 1 (74) 

for I < m, and Mi <m (rf) = for I > m represents the effect of nonperfect detection. 
Since detection of k coincident events can result from various numbers m of 
photons, the conditional state is in general a statistical mixture. Therefore in 
place of equation (|T5D we now have 

Qont(n,k) = J2P<S> P N, V ( n ' m \ k ) \^n,m)(^n,m\, (75) 
m,<E> 

where \ty n , m ) is given in equation fll6|), and P N ^(n,m\k) is the probability of m 
photons being present under the condition that k coincidences are recorded. The 
conditional probability P^^{n,m\k) can be obtained using the Bayes rule as 

P N)V {n,m\k) = -~ — - — — P N ^(k\m) P(n,m). (76) 
P NtV {n,k) 

Here P(n,m) is the prior probability (|T2|) of m photons being present, and ac- 
cordingly, Pjv,7,(n, k) is the prior probability of recording k coincident events, 

PN,«{n, k) = £ P N)V (k\m) P{n, m). (77) 

m 

Second, preparation of the reference mode in a Fock state is a nontrivial 
problem (for a review, see 0; for single-photon Fock states, see also [0, ^5 
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66[; for multiphoton Fock states, see also |)7j |68|, |69|, [7(], fflj). In particular, a 



method for synthesizing multiphoton Fock states from single-photon Fock states 



(produced, e.g., in parametric down conversion) has been proposed |70|) fflf. In the 
scheme, modes prepared in single-photon Fock states are fed into the input ports 
of an array of beam splitters and detectors survey all but one output port so that 
the mode in the free output port is prepared in the sought photon-number state. 
In practice however, it may be more realistic to consider statistical mixtures of 
photon-number states rather than pure Fock states. Let us return to equation 
(|J) and assume that 

Qin = £?inl <8> ftn2, (78) 

where 

ftn2 = 5^Pn \n)(n\ . (79) 

n 

To be more specific, let us consider (as an example of a sub-Poissonian distribu- 
tion) a binomial probability distribution, 

Pn = ("°)p n (l " P) n °~ n if * < n (80) 

and p n = elsewhere (0 <p < 1). Note that for p — > 0, n — > oo, and pno finite the 
binomial distribution QgQP reduces to a Poisson distribution, with n = pno being 
the mean photon number. Using equations (|78|) and ( ff9|) , from equation (|75| ) we 
easily find that after recording k coincident events the conditional mixed state 
now reads 

Qouti(k) =J2p n Qouti{n, k) 

n 

= PnP® PN,r,(n,m\k)ty nym )($ n , m \- (81) 

n,m,<l? 

Accordingly, the probability of detecting the state is the average of Pisr jr] (n, k) 
given in equation (|77D , i.e., 

PnM = Y,PnPN, v (n,k) (82) 

n 

= Yp n P N ^(k\m)P(n,m). (83) 

n,m 

The quadrature- component distributions and the Wigner function of a mixed 
conditional output state ([SI]) are plotted in figure | [P/v i?7 (/c) ~21.4%]. We see 
that the quantum interference is still preserved for realistic values of the number 
of photodiodes and efficiencies (k = 4, N = 20, rj = 90%) and for a sub-Poissonian 
statistics of the input Fock-state mixture (|79|) [n = 3.8, (An) 2 = 0.19]. 
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8 Summary and conclusions 



We have extended previous work on quantum-state preparation via conditional 
output measurement on a beam splitter and shown that when a mode prepared 
in a state |<&) is mixed with a mode prepared in an n-photon Fock state and m 
photons are detected in one of the output channels of the beam splitter, then the 
mode in the other output channel is prepared in a photon-subtracted (n < m) or 
a photon-added (n > m) Jacobi polynomial state which is obtained by applying 
an operator- valued Jacobi polynomial to the state \%) ~ T n |<3>). Since for typical 
classes of input states, such as thermal states, coherent states and squeezed states, 
the states |$) and |\&) belong to the same class of states, the Jacobi polynomial 
states derived from |$) and also belong to the same class of states. Jacobi 
polynomial states are nonclassical states in general, so that subtracting photons 
from or adding photons to them again yields nonclassical states in general. 

The analysis has shown that PSJP and PAJP coherent states can be regarded 
as extremely non-Gaussian squeezed states. A characteristic feature of PSJP and 
PAJP squeezed vacuum states are the well pronounced quantum interferences 
associated with the quadrature-component noise reduction. Moreover, these two 
classes of states represent Schrodinger-cat-like states. Photon-added Jacobi poly- 
nomial states are more nonclassical than photon-subtracted Jacobi polynomial 
states in general. In particular when n = (m > 0) or m = (n > 0), respec- 
tively, then the states reduce to ordinary photon-subtracted or photon-added 
states. Whereas photon-added coherent states are non-Gaussian squeezed states, 
subtracting photons from a coherent state obviously leaves the state unchanged. 

The analysis has further shown that there are close relations to other nonclas- 
sical states that have widely been studied. Hence combining state preparation 
via conditional output measurement on a beam splitter with other schemes of- 
fers novel possibilities of nonclassical-state generation and manipulation, such as 
the generation of near-photon-number eigenstates and squeezed-state excitations. 
Since near-photon-number eigenstates are coherently displaced photon-added co- 
herent states, they can be generated by combining the scheme for photon adding 
with a scheme for coherently displacing a state. The latter can be realized by us- 
ing a second beam splitter whose transmittance is close to unity and which mixes 
the mode prepared in a photon-added coherent state with a mode prepared in 
a strong coherent state. Similarly, a squeezed-state excitation can be prepared 
by appropriately squeezing a photon-added squeezed vacuum state followed by a 
coherent displacement of the state. 

In order to demonstrate the feasibility of generating PSJP and PAJP states, 
we have calculated the corresponding event probabilities. Further, we have also 
allowed for both nonprecise input Fock-state preparation and nonprecise output 
photon counting. For this purpose we have considered sub-Poissonian mixtures 
of Fock states in place of pure Fock states and assumed that photon-chopping 
is adopted for photon counting. The results show that, apart from some smear- 
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ing, typical properties of the states can still be observed even under realistic 
experimental conditions. 
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Appendix A: Proof of equation ([21]) 

To prove equation ([H]), let us consider the operator function 
Ff(n,M; T)=T*X; 

k=0 

which reads in the Fock basis as 



(- 


R 


2 ) k k\ 


(k+ 


v\)\ 



(a.i; 



oo jin 



Ff(n>l;T) = £ 



^o (n+^V- 



i 

fc=0 



r(n+ n+k- 


H)(- 


R 


2\k 


T(k + 


V 


+ 1) 



\n)(n\. 



(A.2) 



Introducing the integral representation of the gamma function T(n + fx + k + 1), 
we have 

F?(n,\v\;T) = Y,-< -/ dtt n ^e 



^0 (n+fi)\ Jo 

' (-\R\ 2 t) k (i s 
j^r(*+M+i) \K 



\n){n\. 



(A.3) 



We now use the sum rule (|2~8D and the integral identity pi 

POO 

/ dtt^e-^Jct) 
Jo 

= r(a)Pf Q - /3 -"- 1) (l - 2c) (Re a > 0) 
and obtain (\R\ 2 = 1 - \T\ 2 ) 



(A.4) 



Ff(n>|;T) 
Z! 



(i+H)iS P ' 



(|^,n+Ai-|H-0 



(2|T| 2 -l)T»H 



Z! 



y(\v\,fl+H-\u\-l) 



2\T\ 2 -l)T n . 



(i + M)i-' v-i-i-v < A ' 5 > 

Combining equations ([18|), fl2~0|), ( |A.1|) and ( |A.5|) eventually yields equation (|2T|). 
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Appendix B: Photon statistics of PSJP and PAJP 

coherent states 



From equation fl2~7|) the photon-number distribution p n ,m(0 = I ( l\^n,m) | 2 °f PSJP 
and PAJP coherent states can be given by 



Pn,m (0 



-^l 2 /!|/3'| 



2/ 



2u 



\R 



2\k 



n 



E 

k=fi 



w 

0(1-3) 



(b.i; 



where 9(n) = 1 for n > and 0(n) = elsewhere. Further, from equations (f24|) 
and (p5|) it can be shown that the antinormally ordered moments of the photon 
number can by given by 



(*n ) m|a P (a t ) P |*»,m> 



l\—2v n n 



ML 



EE(-I^l 2 f +J 



"n,m k=fij=fj. 



(B.2) 



can then be used to 

derive closed solutions for the normally ordered moments of the photon number. 
In particular, writing 



with Xik{ a ) being defined in equation fl32|). Equation 



(n) = (* n ,m|a t a| 1 S r n,m) = (^n,m|aa t |^n,m) 



(B.3) 



and applying equation (|B.2|) , the mean number of photons is calculated to be 

3/1— 2v n n 



(n) 



|/? | " x - , > ; / II \ I II 

k j /u\ /„-\ fti*\V ^X x ) 

x EE 



i=fi i'=fj, 



EE(-I^l 2 ) fc+ ^ 

1)\V) {l - v))\(V - V )\ 



i. 



(B.4) 



In a similar way, closed solutions can also be found for higher-order moments. 
For example, in order to determine the Mandel factor Q = ((n 2 ) — (n) 2 )/ (n) —1, 
knowledge of (n 2 ) is required. It can be obtained by introducing the antinormally 
ordered form 



(n 2 ) = (V n , m \{tfa)*\ib ntm ) 

= (ty la 2 fa tN ) 2 l\l> )— 3(iff laatl^r \ + i 



(B.5) 



and then applying equation (p.2j ). 
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Appendix C: Derivation of equation ( |S4| ) 

According to equation (0), the probability of producing PSJP and PAJP coher- 
ent states is given by 

_ m2 \R\- 2 »n\ 
F{n,m) = e m 



n n 



where 



\T\ 2m m\ 

k — vj \j—u. 



x EEH«l ! )* + 'l™ , m VQW,"), (ci) 



*8&M = £M(™1& (C2) 



k J \ j J (p+v)l 



which for v > (i.e., fi = v and 5 = 0) can be rewritten as 



x 



^l^l 2 )^' 1 liiSrK^+i) )' ( °' 3) 



The term in the curly brackets is nothing but the series expansion of the confluent 
hypergeometric function. We therefore have 

*S (/3 ' ") = h.d$w^ T ^ 3 +1 ' v+1 ' l(3n (c,4) 



Using the relations [61 



$(a, b, z) = e z $(b - a, b, -z), 



we then obtain 

)k 



*!»> = V wrt Prwui-vr?) 



j! ^/ « Wl 2 )' 

The form ( |35D of Xk,jWi u ) use d m equation ([34]) follows by applying standard 
formulas for derivatives of Laguerre polynomials |H[ . 
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For v < (i.e., \x = and 5 = —v) equation ( |C.2|) can be rewritten as 



n=0 



k I \ j I n\ 



1 dk m w v ( n - v+ A H^L (c 7) 



Using the relations |61 



Li (°) 

\ J J 



n=0 U - 



we find that 

%^^^r^ r( -^y (c.9) 

By again applying standard formulas for derivatives of the Laguerre polynomials 
|6lH equation ( |C.9| ) takes the form used in equation (|35| ) . 
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BS 




\n) 

Figure 1: Scheme of the experimental setup. When a signal mode prepared in 
a state £ inl is mixed (beam splitter BS) with another input mode prepared in a 
Fock state |n) and in one of the output channels of the beam splitter m photons 
are recorded (detector D), then the quantum state g out i of the mode in the other 
output channel "collapses" to either a PS JP state (n <m) or a PAJP state (n > 
m). 
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Figure 2: The probability of producing (n — 2, m = 3) PSJP coherent states 
(dashed line) and (n = 3, m — 2) PAJP coherent states (solid line) is shown as a 
function of \(3\ for two values of the beam-splitter transmittance [(a) |T| 2 = 0.4; 
(b) |T| 2 = 0.81]. 
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Figure 3: The quadrature-component distributions of (a) PSJP and (b) PAJP 
coherent states for (3' = 2.07 (|/3|=2.3, |T| 2 = 0.81) [(a) n = 2, m = 3, (b) n = 3, m 
= 2]. 
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Figure 5: The Wigner function of (a) PSJP and (b) PAJP coherent states for (3' 
= 2.07 (|/3|=2.3, |T| 2 = 0.81) [(a) n = 2, m = 3, (b) n = 3, m = 2]. 
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Figure 6: Experimental setup for preparing near-photon-number states (crescent 
states). A mode prepared in a coherent state \a) is first mixed (beam splitter BSl) 
with a mode prepared in a Fock state \n) and a zero-photon output measurement 
is performed (detector D). The output mode prepared in a photon-added coherent 
state {cv) n \2f3) (2j3 = Ta) is then mixed (beam splitter BS2 with transmittance 
T close to unity) with a mode prepared in a strong coherent state \ciil) such that 
P = -(R/f)a L . 
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Figure 7: Experimental setup for preparing squeezed-state excitations. A mode 
prepared in a squeezed vacuum state S(() |0) is first mixed (beam splitter BSf) 
with a mode prepared in a Fock state \n) and a zero-photon output measurement 
is performed (detector D), so that the output mode is prepared in a photon-added 
squeezed vacuum state (a^) n 5 , (2^|0) [k(2£) =T 2 k(()}, which is then squeezed, e.g., 
by a degenerate parametric amplifier (DPA) to realize the squeezing parameter 
— £. Finally the so squeezed output mode is mixed (beam splitter BS2 with 
transmittance T close to unity) with a mode prepared in a strong coherent state 
|a£,) such that (3 = (R/T)ctL. 
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Figure 8: The quadrature-component distribution (a) and the Wigner function 
(b) of a mixed conditional output state (^) for a coherent input state with j3' 
= 2.07 (\/3\ =2.3, |T| 2 = 0.81). The parameters of the photon-chopping detection 
scheme are k = 4, N = 20, rj = 90% and those of the input Fock-state mixture (fT9|) 
are Uq = 4 and p = 0.95. 
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